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DISCOVERY OF A PULSAR WIND NEBULA CANDIDATE IN THE CYGNUS LOOP 
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ABSTRACT 

We report on a discovery of a diffuse nebula containing a pointlike source in the southern blowout 
region of the Cygnus Loop supernova remnant, based on Suzaku and XMM-Newton observations. The 
X-ray spectra from the nebula and the pointlike source are well represented by an absorbed power-law 
model with photon indices of 2.2±0.1 and 1.6±0.2, respectively. The photon indices as well as the flux 
ratio of -F n ebuia/-Fpointiike ~ 4 lead us to propose that the system is a pulsar wind nebula, although 
pulsations have not yet been detected. If we attribute its origin to the Cygnus Loop supernova, then 
the 0.5-8 keV luminosity of the nebula is computed to be 2.1xl0 31 (ci/540pc) 2 ergss -1 , where d is the 
distance to the Loop. This implies a spin-down loss-energy E ~ 2.6 x 10 35 (d/540 pc) 2 ergs s -1 . The 
location of the neutron star candidate, ^2° away from the geometric center of the Loop, implies a high 
transverse velocity of ^1850 (9/2°) (d/540pc) (i/10kyr) -1 kms -1 , assuming the currently accepted 
age of the Cygnus Loop. 

Subject headings: ISM: individual objects (Cygnus Loop) — ISM: supernova remnants — pulsars: 
general — stars: neutron — stars: winds, outflows — X-rays: ISM 



1. INTRODUCTION 

The Cygnus Loop supernova remnant (SNR), the X- 
ray image of which is shown in Fig[TJ is one of the 
brightest SNRs in the X- ray sky. Because of its prox- 
imity, 540(+100,-80)pc (jBlair et all 120091) . it is an ex- 
tremely important object that is often called a proto - 
typical middle-aged SNR (lOkvr: ILevenson et al.lll998l ). 
The X-ray morphology is an almost perfect circular shell 
except for a southern blowout region. The origin of the 
blowout has been a matter of debate; it may be caused by 
either low ambient density (Aschenbach & Leahy Il999t 
Uchida et all 120081) o r a second SNR (jUvaniker et al l 



2002: ISun et a l. 2006). 



Another mystery for the Cygnus Loop is the absence 
of a central compact remnant. It is believed that the 
Cygnus Loop is the result of a core-collapse SN, be- 
cause the blast wave is now hitting the walls of the 
cavity that was most likely create d by a strong stel- 
lar wind from the progenitor (e.g.. [Charles et al l 119851 : 
IHester et all 119941: ILevenson et al.1 119971) . The compar- 
atively small size of the cavity (R ~ 13 p c at a dis- 
tance of 540 pc) led ILevenson et al.1 (|1998f ) to suggest 
that the progenitor star was of spectral type later than 
B0, ^15M . This view is further supported by recent 
X-ray abundance measurements which h ave led to pro- 
genitor mass estimates of 12-15M(T) (e.g. . iTsunemi et al.l 
2007; Ki mura et al.l l200l lUchida et all [20111) • Such a 
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progenitor star should have formed a neutron star during 
the SN explosion. Although considerable effort has been 
devoted to searching for a neutron star in the Cygnus 
Loop over nearly three decades , none has been found yet 
(e.g.. IMivatalt^lT99l 12001) . 

Here, we report on the discovery of a possible pul- 
sar wind nebula (PWN) in the southern blowout region 
of the Cygnus Loop, based on X-ray observations with 
Suzaku and XMM-Newton. This Letter focuses on the 
discussion about the PWN candidate, while other results 
based on these data have been published elsewhere (e.g., 
Tsunemi et all 2007; U chida et all 120111 : [Katsudaet_a]J 
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2. OBSERVATIONS 

We have conducted over 80 pointing observations of the 
Cygnus Lo op using the Suzaku X -ray Imaging Spectrom- 
eter (XIS: iKovama et"aH I2007D and the XMM-Newton 
European Photon Imagin g Camera (EPIC: [Turner et al.l 
l2001t IStriider et "aT1l2001t ) , covering nearly the entirety of 
this large SNR. The fields of view (FOV) of these obser- 
vations are overlaid on the R OS AT all-sky survey image 
in Fig.[TJ We focus here on one XIS/EPIC observation lo- 
cated in the southern blowout region. The XIS and EPIC 
observations were performed on 2011-05-07 (Obs.ID: 
506013010) and 2006-05-15 (Obs.ID: 0405490301), re- 
spectively. 

For the Suzaku XIS data, we use cleaned event data 
prepared by the Suzaku operations team. The net expo- 
sure time is 60.3 ks for each XIS. The XMM-Newton data 
suffer severely from high background (BG) flares due to 
soft protons throughout the observation. Nonetheless, 
the EPIC image may provide us with useful information 
on spatial structures, thanks to much better angular res- 
olution than the XIS. We thus use relatively clean time 
regions, where the count rates in the 5-12 keV are less 
than 5ctss _1 for MOS1/2 or 50ctss _1 for pn. The ef- 
fective exposure times obtained are 7.8 ks, 8.0 ks, and 
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Figure 1. ROSAT all-sky survey image of the Cygnus Loop, 
scaled as the square root of the surface brightness. The energy 
band used is 0.2-2 keV. All of the Suzaku XIS and the XMM- 
Newton EPIC fields observed to date are overlaid as boxes and 
circles, respectively. The XIS and EPIC images of the yellow FOV 
is shown in Fig. [2] The geometric center of the Cygnus Loop is 
indicated by a cross. 

6.0 ks for MOS1, MOS2, and pn, respectively. All the raw 
data are processed using version 11.0.0 of the XMM Sci- 
ence Analysis Software. Further analyses are done with 
heasoft tools of version 6.11.1 and the latest CALDB 
files updated on 2011-11-09. 

3. ANALYSIS AND RESULTS 

Through our comprehensive analyses of the X-ray data, 
we have discovered a hard X-ray-emitting nebula in the 
southern blowout region of the Cygnus Loop, i.e., the 
XIS/EPIC FOV illustrated in yellow in Fig. [TJ Figure H] 
left and center are vignetting-corrected XIS images in 
0.5-1 keV and 1-10 keV, respectively, for which non X- 
ray BG is subtrac ted by using the xisnxbgen software 
(|Tawa et al.l 12008ft . While the left panel is dominated 
by thermal emission from the Cygnus Loop, the central 
panel represents mixture of cosmic X-ray BG and hard 
X-rays from astrophysical sources. In the hard-band im- 
age, we see a bimodal diffuse feature — our target of in- 
terest in this Letter. An XMM-Newton image in Fig. [5] 
right successfully resolves it into a northern diffuse neb- 
ula and a southern pointlike source. It should be noted 
that the short exposure time and the presence of soft 
proton events in the XMM-Newton image make it diffi- 
cult to detect faint diffuse emission between the pointlike 
source and the northern nebula. 

The position of the pointlike source is determined to be 
[RA, Dec] = [20:49:20.309, +29:01:05.57 (J2000)], by us- 
ing the emldetect software. The statistical uncertainty 



is negligible compared with the astrometric uncertainty 
of 2" (b ased on the XMM -Newton Calibration Technical 
Note - iGuainazzi 1 120111) . The source has been identi- 
fied as 2XMM J204920.2+290 106 in the XMM-N ewton 
Serendipitous Source Catalog (Watson et al. 2009j). No 
obvious optical and infrared counterpart is found in 
the Palomar Observatory Sky Survey and 2MASS All 
Sky Survey, respectively, although there is an appar- 
ent compact galaxy whose peak is ~3" away from the 
X-ray pointlike source and its optical emission extends 
to the source position. The R-band upper limit is es- 
timated to be R>19 mag. No radio counterpart is de- 
tected in an archival VLA 1.4 GHz image. The radial 
profile of the pointlike source is found to be consistent 
with the EPIC point-spread function at the same off- 
axis angle. The northern nebula is positionally consis- 
tent with an infrared object 2MASX J20491447+2903237 
which is identified as an extended extra galactic source 
(jSkrutskie et al.ll2006l ). while its extent appears to be 
quite small and thus there is no evidence for the physical 
association to the X-ray nebula. 

We examine the XIS spectra of the pointlike source and 
the surrounding nebula. The spectral extraction regions, 
a l'-radius circle around the pointlike source and a poly- 
gon tracing the edge of the diffuse nebula, are shown 
in Fig. [5] center. BG is taken from an elliptical annu- 
lus shown in Fig. [2] left and center. The BG produced 
in this way includes not only usual X-ray BG but also 
thermal emission from the Cygnus Loop itself (typically, 
kT c -0.3 keV and n c t ~10 n cm~ 3 s: lUchida et alll200l . 

The XISO+3 spectra together with the area-normalized 
local-BG spectra are shown in the upper panels of Fig. [3j 
and the local-BG subtracted spectra are plotted in the 
lower panels. The photon numbers after BG subtraction 
are 753 and 4534 for the pointlike source and the nebula, 
respectively. Since emission below 1.5 keV is dominated 
by Cygnus Loop's thermal emission, it is in principle dif- 
ficult to subtract the local BG properly in the soft X- 
ray band. We thus use the 1.5-10keV band for spectral 
analysis. In this case, the intervening column density 
cannot be constrained, so we tentatively fix the value of 
iVn t o 4x 10 20 cm" 2 wh ich is typical for the Cygnus Loop 
(e.g.. Hnoue et al.lll980f ). Before fitting, each spectrum is 
grouped into bins with at least ~50 counts in prior to 
background subtraction, which allows us to perform a 
X 2 test. In the calculation of \ 2 values, we use the stan- 
dard weighting method, i.e., the square root of detected 
counts. 

We find that either a power-law m odel or an optically 
thin thermal emission model (apec: iSmith et all .2001) 
can adequately describe the data. The thermal emission 
model, however, requires a high electron temperature of 
2-4 keV, and very low abundances of <0.8 times the solar 
value. Such a high temperatu re has never been rep orted 
from the Cygnus Loop (e.g., iTsunemi et al.l 120071 ) . and 
the low abundances are unusual in astrophysical sources. 
We thus reject the thermal model from an astrophysical 
point of view. Also, a blackbody model is safely rejected 
from a statistical point of view. Adding a blackbody 
or thermal emission component to the power-law com- 
ponent does not improve the fits, either. In this way, 
we conclude that the power-law is the most reasonable 
model to describe the data. The best-fit parameters from 
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1 Suzaku XISO+1+3 ( 1 - 1 Oke V) I -XMM-Hewton 
I _ f , ,. t /\ MOSl+2(0.2-10keV) 




2XMM J204920.2+290106 



Figure 2. Left: Soft-band (0.5— IkeV) XIS1 image of the yellow box in Fig. [T] The image is scaled as a square root of the surface 
brightness and corrected for vignetting effects after subtraction of non X-ray BG. A white elliptical annulus shows where we extract a local 
BG. Center: Same as left but in the hard band. The data taken by XISO, XIS1, and XIS3 are summed to improve photon statistics. The 
spectral extraction regions are shown as a white circle (for the pointlike source) and a polygon (for the nebula). Right: Closeup image 
of the white box (6' square) in Fig. [2] center taken by the XMM-Newton MOS1+2. BG is not subtracted and vignetting effects are not 
corrected. 



0.1 
0.01 

> 

X, io- 3 



10-- 



The pointlike source 



i #]+++++++ 

-Witt +++ Source (XIS0+3) 

tft , ++_r_+ + 



Local BG (XIS0+3) 



3) 



o.i 



0.5 



10 



0.01 



U 1Q-* 



Energy (keV) 

The pointlike source 



XIS0+3 
XIS1 




2 5 
Energy (keV) 



Jl.. , The nebula 



0.01 



W3 



fttttt, Source (XIS0+3) 



Local BG (XIS0+3) 



io-' 



0.5 1 2 5 

Energy (keV) 

The nebula 

tf%j£bki. XIS0+3 
. - XISl 




Energy (keV) 



Figure 3. Left upper: Suzaku XIS (0+3) spectrum for the pointlike source (black) and for the area-normalized BG spectrum (gray). 
Right upper: Same as left but for the diffuse nebula. Left lower: Local-BG subtracted XIS spectra for the pointlike source along with the 
best-fit power-law model. Black and red correspond to the FI (XISO+3) and the BI (XISl), respectively. Lower panel shows residuals. 
Right lower: Same as left but for the diffuse nebula. 



the power-law model are summarized in TableQ] We note 
that consistent results are obtained if we let the soft X- 
ray band (i.e., 0.5-1.5 keV) remain and fit the spectra by 
a power-law plus local-BG model whose normalization is 
allowed to vary freely. Also, no significant spatial varia- 
tion of the photon index in the synchrotron nebula has 
been detected owing to insufficient photon statistics. 

To search for long-term spectral variability of the 
pointlike source, we also fit the EPIC spectra. While the 



Table 1 

Power-law spectral-fit parameters 



Region 


r 


Flux" 


X 2 /cI.o.f. 


The pointlike source 


1.6±0.2 


2.0±0.5 


6.6/12 


The nebula 


2.2±0.1 
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96.1/86 


Note. — Errors quoted a 
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a TJnabsorbed 0.5-8 keV flux 
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data are affected by soft protons, flux estimates would be 
relatively robust if we subtract a local BG. We find that 
the fluxes (and photon indices) inferred from the MOS 
and pn spectra are marginally consistent with those in- 
ferred using the XIS. In addition, a Fourier search does 
not give significant pulsations. This might be due to 
poor photon statistics as well as insufficient time res- 
olution (0.073 s, 2.6 s, and 8 s for pn, MOS, and XIS, 
respectively) . 

4. DISCUSSION 

We have presented the discovery of a diffuse neb- 
ula in the Cygnus Loop containing a pointlike source. 
The spectra of both are well represented by a power- 
law model, reminding us of a PWN. A PWN ori- 
gin is further suppo rted by the follow i ng o bser- 
vational features. Kargalts ev fc Pavlovl (|2008l ) re- 
ported a strong correlation between PWN luminosi- 
ties and the nonthermal luminosities of pulsars, namely 
£pwn (0.5-8 keV)/£pSR (0.5-8 keV) is nearly constant at 
4. Our measured £ nebula (0.5-8 kcV)/^ P ointiikc (0.5-8 kcV) 
of ^3.8 is in good agreement with this correlation. 
Above, we assume that the nebula extends over the 
pulsar candidate with its mean surface brightness, re- 
sulting in Nebula (0.5-8 kcV) ~ 5.7x 10 -13 ergs cm" 2 s -1 
and i^pointiikc (0.5-8 kcV) ~ 1.5 xlO -13 ergs cm -2 s -1 ._ Fur- 
thermore, the measured photon indices are consistent 
with those of PWNe, and the nebula's photon index is 
too hard to be interpreted as shell emission from rel- 
ativistic particles accelerated at an SNR shock (gen- 
erally, T ~ 2.5 — 3). It should be also noted that 
the spectrum of the nebula is softer than the pointlikc 
source, typical of PWNe. From a morphological point 
of view, the brightest part of the nebula is offset from 
the pointlike source. This is unusual for PWNe, but 
complicated structures are often found in PWNe (e.g., 
iKargaltsev fc Pavlovl 120081 ): hence the morphology can- 
not be a strong reason to eliminate the PWN hypoth- 
esis. Also, the pointlike source is not likely to be a 
background active galactic nucleus or a foreground nor- 
mal star, since it has no obvious radio or optical coun- 
terpart and its X-ray spectrum is too hard for a nor- 
mal star. While low mass X-ray binaries (LMXBs) have 
high X-ray to optical luminosity ratios, this possibility is 
ruled out from the low X-ray flux of the pointlike source, 
^-lokov ~1.3xl0 -13 ergscm -2 s -1 based on the XIS; 
a normal LMXB luminosity of ^-TO 36 ergss -1 requires a 
distance of 250 kpc which is far outside the Galaxy. On 
the other hand, these properties are consistent with the 
expectations for an isolated neutron star, even though 
the conclusive evidence of a pulse period is lacking. Nev- 
ertheless, we propose that this system is a PWN 

We note, however, another possibility: that the sys- 
tem consists of a galaxy cluster and an unknown point- 
like source which might not be related to the clus- 
ter. A brief justification can be made by comparing 
the nebula's properties with those of galaxy clusters. 
A Lx-temperature c orrelation of galaxy clusters (e.g., 
iFukazawa et al.|[2004f) would result in nebula's luminos- 
ity, L 2 -iokeV = 2.0 x 10 44 (fcT/3keV) 2 79 ergss -1 . Com- 
bining the luminosity with nebula's flux, i^-iokoV = 
2.7 x 10 -13 ergss -1 , we obtain a distance of ^2.5 Gpc. 
The apparent radius of 3' can be then translated to a 



real radius of ^2 (d/2.5 Gpc) Mp c, which is typical fo r 
galaxy clusters with kT of 3 keV ( Fukaza wa et al.ll2004l) . 
Therefore, the X-ray luminosity, temperature, and size 
are all in agreement with those of galaxy clusters. In 
addition, based on the ROSA T PSPC survery of galaxy 
clusters (jBurenin et al.l 120071 ). a surface density of clus- 
ters is ^0.3 deg -2 at ^2.8xl0 -13 ergs cm -2 s -1 which is 
the nebula's brightness in 0.5-2 keV. Given the Cygnus 
Loop's large dimension of ^10 degree 2 , we expect ~3 
clusters inside the Loop. Thus, a chance probability to 
detect a galaxy cluster in the Loop is fairly large. 

While other possibilities have not been excluded, this 
system is most likely a PWN based on its observed 
properties. We thus further discuss the PWN scenario. 
Whether it is physically associated with the Cygnus Loop 
is an important point. Evidence favoring a common 
origin is the fact that, even though the Cygnus Loop 
should have originated from a core-collapse SN, our sur- 
vey of the Suzaku and XMM-Newton data has not de- 
tected any PWN or neutron star except for the PWN 
candidate disclosed here. If the objects are associated, 
then the distance to the PWN s hould be the sam e as 
the Cygnus Loop, i.e., d ~540pc (jBlair et al.ll2005T) . If 
this is the case, the 0.5-8 keV luminosities of the pulsar 
candidate and the nebula are computed to be 5xl0 30 
(d/540 pc) 2 ergss- 1 and 2.1xI0 31 (d/540pc) 2 ergss -1 , 
respectively. This would be one of the lowest lumi- 
nosities measured from the ~50 known X-ray PWNe 
(Kargaltsev & Pavlov 2008). The physical size of the 
nebula would be ^0.5 (# r /3')(d/540pc) pc, where T is 
the angular radius of the nebula. This size is an order 
of magnitude smaller than e xpected based on the PWN 
size-age relation reported bv iBamba et al.l (|2010f ) . This 
discrepancy is mitigated if the PWN is farther away than 
the Cygnus Loop. Therefore, the PWN candidate is ei- 
ther a distant, ordinary PWN unrelated to the Cygnus 
Loop or a nearby, peculiar PWN. This point should be 
revisited by future detailed X-ray spectroscopy that will 
allow for an Nn measurement. 

We next estimate some important par ameters, as is 
usual l y performed for putative PWNe Ce.g-.lHughes et al.l 
[20011 iQlbert et al.ll200l IGaensler et al.ll2003D . A well- 
known relation between PWN luminosity and spin- 
down loss power, 77 = Lx/E ~ 8 x 10 -5 (most re- 
cently, iVink et all 1201 ID . indicates that E ~ 2.6 x 
10 35 (ci/540pc) 2 ergss -1 . The estimated E allows us to 
deduce a pulsar period of P ~ 0.48 [2/{n- l)] ' 5 (E/2.6 x 
10 35 ergss -1 ) -0,5 (i/10 kyr) -0 5 s, assuming a standard 
moment of inertia (/ = 10 45 gcm 2 , i.e., a standard 
neutron star radius of 10km and a mass of 1.4 M Q ), 
a braking index of 3, and a negligible initial pe- 
riod. Also, the period derivative can be estimated 
to be P = P/[(n - l)t] - 7.7 x 10 -13 (25/2.6 x 
10 35 ergss -1 ) -0 ' 5 (i/10kyr) -1 ' 5 ss -1 , leading to a sur- 
face magnetic field estimate of B ~ 1.9 x 10 13 (i?/3.4 x 
10 35 ergss -1 ) -a5 (i/10kyr) -1 G. Comparing these pa- 
rameters with those of the other X-ray-emitting PWNe 
(jKargaltsev fc Pavlovl 12001 . we find that E is some- 
what smaller than others derived at the age of lOkyr 
and that the P and B values are among the largest of 
all the listed PWNe. Alternatively, noting that there is 
a correlation between E and r, E = 41.1 — 1.08 log r 
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for r < 10 4 kyr, we can infer r for the estimated E of 
2.6 xlO 35 (d/540pc) 2 ergss~ 1 to be -180 kyr. This would 
lead to P ~ 112 ms a nd B ~ 10 12 G, which are com- 
mon for X-ray PWNe (|Kargaltsev fc Pavlov! [2008h . We 
also point out that the estimated E implies a gamma- 
ray luminosity for a possible pulsar of 1.6x1 (E/2.6 x 
10 35 e rgs s" 1 )" - 5 ergs s" 1 in 0.1-100 GeV l|Abdo et al.l 
120101 ). This should be easi ly detected with the Fermi 
LAT, which is not the case (Katagiri ct al. 20"TTI). 

It is interesting to note that the neutron star 
candidate is located far from Cyg nus Loop ' s geo - 
metric center defined by Einstein (|Ku et al.l 119841) . 
as can be seen in Fig. [TJ If we assume that 
the pointlike source started moving from the ge- 
ometric center lOkyr ago, its transverse velocity 
is -1850(6>/2°)(d/540pc)(t/10kyr)- 1 kms- 1 , and its 
proper motion ~0".72 (6»/2°)(t/10kyr)" 1 yr" 1 , where 
is the angular distance between the neutron star candi- 
date and the geometric center and t is the best-estimated 
age of the Cygnus Loop. Such a high velocity is one of th e 
fastest of known neutron stars (e.g., iKaspi et al.l 120061) . 
However, we should keep in mind that the velocity esti- 
mate is quite uncertain. The explosion location is quite 
uncertain (and thus 9), because the SN went off in a cav- 
ity In fact, it has been propo sed that the SN occurred i n 
the southern blowout region (|Tenorio-Tagle et al.lll985l) . 
The age is also a matter of debate; it might not be accu- 
rate to 50%. 

We have checked that the existing data including Ein- 
stein, ROSAT, and ASCA do not allow us to measure 
the proper motion of the source; the source is not de- 
tected, possibly due to insufficient exposure times as well 
as low detection efficiency in the hard X-ray band. It is 
also difficult to measure proper motions using the XIS 
data, given the larg e position uncertainty of ~1' due to 
thermal fluctuation (jUchivama et al.H2008|) . Thus, high- 
resolution images with Chandra or XMM-Newton will be 
required. Confirmation of a large proper motion towards 
the south would be direct evidence that the pointlike 
source is indeed the neutron star remnant of the Cygnus 
Loop. In this case, we expect a com etary bowshock struc- 
ture around the neutron star (cf. Gaensler c t alJ[2004T ) 
due to the high Mach number of ~6, given a reason- 
able sound speed of c s = (jkT) / (/xm p ) ~ 300 km s -1 , 
where 7 is a specific heat ratio (5/3), kT is the plasma 
temperature (0.3 keV), and \x is the mean molecular mass 
(~0.6 for solar abundance plasmas), and m p is the pro- 
ton mass. Note that such a bow shock would be directed 
away from the geometric center of the Loop and thus 
not necessarily related to the northern nebula. A much 
higher Mach number would be expected if the system is 
already outside the hot plasma and is now proceeding 
into cold interstellar medium. On the other hand, if the 
proper motion turns out to be small, the location of the 
PWN candidate near the center of the blowout appears 
to support the suggestion that the s outhern blowout is a 
separate SNR (|Uvaniker et al.ll2002D . 

5. CONCLUSION 

Using Suzaku and XMM-Newton, we discovered a 
hard X-ray-emitting diffuse nebula containing a point- 
like source in the southern blowout region of the Cygnus 
Loop. Their properties suggest that it is most likely a 



PWN. The physical relation to the Cygnus Loop is still 
uncertain at this point. Future observations are essential 
in order to confirm that the object is a PWN as well as 
to clarify whether it is associated with the Cygnus Loop. 
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